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I. INTRODUCTION 
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The study of influence of initial inhomogeneities upon the evolution of cosmological models is an important issue 
in Cosmology, both in order to understand the formation of large-scale structures as well the smoothing away of 
3t ■ these inhomogeneities. One way to deal with the inhomogeneities is to consider perturbations of known solutions, 
like Friedmann- Robertson- Walker models or homogeneous Bianchi type solutions. Another way is to assume inho- 
mogeneous solutions from the beginning and study their dynamical evolution; it is this that we are concerned with 
here. The simplest inhomogeneous models are the so called diagonal G2 cosmologies for which the spacetime admits 
two commuting Killing vectors whose orbits are 2-dimensional space-like surfaces. They represent inhomogeneities of 
spatially homogeneous models and can be considered as gravitational waves of a single polarization propagating over 
a homogeneous background Q] . Using a few solution generating techniques a large number of exact solutions with 
different sources have been found 0]- 
' The asymptotic evolution of homogeneous Bianchi type models and both their behavior near the initial singularity 
as well as their future state have been widely studied Q. One of the facts relevant to our paper is that at late 
y—i ' times Bianchi models can be described as self-similar solutions. However, the asymptotic behavior of inhomogeneous 
Cn| , solutions, particularly G2 metrics, revealed a task much more difficult than that of homogeneous models. Contrary to 
the homogeneous metrics for which the Einstein equation reduce to an autonomous system of differential equations 
that can be analyzed using techniques from the theory of dynamical systems, the field equations of inhomogeneous 
^-r . spacetimes are partial differential equations. In ref. [3] the evolution of a particular class of G2 metrics is studied 
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showing that in that case all the solutions are asymptotically self-similar. The same result was obtained in £| and in [§J 
where special families of G2 solutions were considered with a scalar field. On the other side, it has been suggested that 
fluctuations might evolve from arbitrary initial conditions to self-similar form [6j . From all these results it is reasonable 
1 , to regard the self-similar solutions as describing the long time asymptotic of inhomogeneous metrics. Hence it would 
£~ ' ■ be worthy to analyze the asymptotic behavior of general class of G2 metrics by using a method that emphasizes the 
rr| scaling properties of the underlying field equations. 

In this paper we study the asymptotic evolution of vacuum G2 metrics using a different approach than that used 
in the above referred papers: we make use of renormalization group (RG) tools to study the structurally stable 



characteristics of the inhomogeneous metrics. Recently, RG techniques have been exhibited as a powerful implement 
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to study the asymptotic behavior of partial differential equations [Tj, |8|, |9|, [10j . The RG method has been applied to 
study homogeneous and isotropic universe, an spherically symmetric dust collapse [Tl| , critical phenomena related with 
gravitational collapse 0], Newtonian cosmology 0, homogeneous flat causal bulk viscous cosmological models Q 
and the theory of perturbations of isotropic universe with dynamically evolving Newton constant and cosmological 
constant [HE3. 

The plan of the paper is the following: in Section^we illustrate the application of the RG method to a homogeneous 
Bianchi type metric with a scalar field. This case has been studied before and we recover the attractors of the system 
by means of the RG method. In Section ITTT1 we apply the RG technique to the vacuum, diagonal, G2 metric. We find 
the fixed-points and analyze their stability as well. We conclude with Section llVl 



II. THE RG METHOD: AN ILLUSTRATIVE EXAMPLE 



It is well known that the asymptotics of partial differential equations can often be found from the consideration of 
scaling solutions (the equivalence of RG theory and the theory of intermediate asymptotics was shown by Goldenfeld 
et. al. 17]). Though it is usual practice to find the similarity variable to analyze scaling solutions from a combination 
of variables using dimensional arguments, however, there is a large class of problems where this cannot be done (l8| . 
RG provides a systematic approach for finding the scaling variables as well as the asymptotics of partial differential 
equations. 



2 




The general (RG) method adapted to partial differential equations is well known 0. Therefore, instead of repro- 
ducing the whole procedure, we explain this method using an illustrative example which has been studied with other 
methods and has a non-trivial asymptotic structure. We consider a class of anisotropic cosmological models given 

by El 

ds 2 = -dt 2 + a(t) 2 dx 2 + b(t) 2 e 2mx dy 2 + c{t) 2 e 2x dz 2 . (2.1) 

This represents a one parameter (to) family of Bianchi models with Bianchi III, V and VIo for, respectively, to = 0, 1 
and —1, and Bianchi Vlh (h = m + 1) for all other vlaues. 

For a direct comparison with earlier works we choose the same variables, namely, shear (a) and expansion (9) for 
our analysis. Note, however, that the feasibility of the earlier analysis depends crucially on finding the "right" set 
of variables to work-with, otherwise system is too complex to analyze completely, even in this homogeneous case. 
Moreover, the choice of these variables is model dependent and rather ad-hoc. Therefore, one also needs to be lucky 
to be able to analyze. Whereas, with RG one can directly work with the metric functions, as will be done in the next 
Section. 

The shear and expansion for this spacetime are given by 

(2.2) 
(2.3) 

Here an over-dot signifies derivative with respect to t. We consider energy momentum tensor describing a minimally 
coupled homogeneous scalar field (j)(t), given by 

= - 9fiU Q<^' Q + F(0)J , (2.4) 

with an exponential potential V(<p) = Acxp(fc0). Here, A > and k are both constants. The evolution equations are 
then given by 

9 = -2a 2 - 9 - - ft + V(cf>) , (2.5) 

a = -a9+ p(m) (V - 3<r 2 - ^<j> 2 - W((j>)j , (2.6) 

4> = -6j>-kV(4>) , (2.7) 
V = k<j)V , (2.8) 

where pirn) = (1 — m) / (3-\/3(l + m + to 2 )). Since above system is an autonomous system of differential equations, in 
the previous study [19j the theory of dynamical systems was used to analyze its asymptotic behavior using "expansion- 
normalized variables". We will use, as an alternative method, the RG technique and recover the same results. Again, 
our goal is to illustrate the method in this section and is applied it to a new scenario in next section. 

We find it useful to work with a compact notation and, therefore, define a new indexed variable Uiit) that signifies 
the set {9, a, <f>, V}, for i = 1,2, 3, 4, respectively. Our interest is in the asymptotics of solution of the form 

lim Ui (t) = r°X(l) , (2.9) 

t — >oo 

where argument "1" signifies initial value of the quantity. It is convenient to fix initial time as t = 1. Now we 
will illustrate how RG method gives a systematic procedure to fix oti and determine the scaling function w*(l) as 
fixed-points of RG equations (in a inhomogeneous case the argument of u* will not be a constant but will depend on 
a combination of time and spatial variable, the scaling variable) . 
Let us consider transformations of the form 

t — ► Lt , 
u^t) — > U^t) = L ai Ui(Lt) , 

(2.10) 
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which leave the equations invariant, i.e., if Ui(t) is a solution so is Ui(t). We use a number L > 1 as a parameter 
of scale transformation. Note that unlike in the application of RG method to quantum field theories or statistical 
mechanics, there is no natural way to choose a scale L here. Moreover, L never appears explicitly in the equations. 
From equation set (|2.5|) - (|2.8[l and (|2.1UI) it is straightforward to recover 

ai = a 2 = a 3 = — = 1 . (2-11) 

Now, defining L = exp(r), and using Eq. I|2.10(l in (|2.5|I - H2.8J) we get a new set of evolution equations: 

— = ai U i + —. (2.12) 

The procedure of defining r is, in a sense, analogous to summing over all degrees of freedom corresponding to 
fluctuations of scale less than L and then re-scaling everything by L _1 . This new set of evolution equations define 
the RG transformations. Since scaled quantities satisfy same evolution equations the time derivative in the equation 
set above can be replaced with equations l|2.5ll - 12.8|l . and we have 
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(2.13) 

3r, ) . (2.14) 

(2.15) 
(2.16) 



All the quantities on the right hand side are evaluated at t = 1. Scale-invariant solutions emerge now from the 
fixed-point structure of the RG map, which is defined by 

^=0, (2.17) 
dr 

with U* being the fixed-points. The complete set of fixed-points for this system are given in the table below : 
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Table I: Fixed-points of Bianchi Vlh models 



Therefore, we have four exact scale invariant solutions of the form l|2.9|l with exponents ai given by l|2.11|l and uUl) 
being the fixed-points given above. Since the stability properties of these solutions is well studied in literature [19j . 
we move now to the more general inhomogeneous case. We analyze the fixed point structure as well as their stability 
properties, which is not done before. 



III. INHOMOGENEOUS CASE 



We begin with a family of metrics which have been a very useful tool for stud y in g inhomogeneous cosmological 
models, known as the generalized Einstein- Rosen spacetimes (see for example, |20|,|2l|). The line element is of the 
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form: 

ds 2 = f(dz 2 - dt 2 )+-/ ab dx a dx b , a, b =1,2. (3.1) 

Here x 1 — x, x 2 — y, and both / and ^ a b are functions only of z and t. Metric l|3.1[l is fairly general and includes 
Bianchi type models I to VII. In this paper we are interested in vacuum solutions only. We can impose now 

det 7ab = t 2 , (3.2) 

since in vacuum case one of the field equation takes the form 

(det lab )% /2) - (det 7ah )^ /2) = . (3.3) 

t and z are two independent solutions of this equation. Under this condition the coordinates t and z are called as the 
canonical coordinates; and there is no loss of generality in this choice |20|. To simplify the analysis we now specialize 
to the diagonal metrics, i.e., metrics with a single polarization. 
The spacetime can be written in the form 

ds 2 = f{dz 2 - dt 2 ) + t(h 2 dx 2 + hT 2 dy 2 ) , (3.4) 

where / and h are functions of t and z only. These metrics admit an Abelian G2 group of isometries with 2 spacelike 
commuting killing vectors d x and d y . The set of Einstein field equations in the vacuum case reduces to 

h « = (3 - 5) 

Eqs. H3.5J) and (|3.6|) are the evolution equations for h and / , respectively, and it is easy to check that the remaining 
field equations are identically satisfied. We now follow the prescription given in the previous section to recover the 
exact scaling solutions. 

Let us consider the following scale transformation 

z — > Lz , 
t — > LH , 
h(t,z) — ► 4>(t, z) — L a h(L e t, Lz) , 

f(t,z) — » i>(t,z) = LPf(L%Lz). (3.7) 

Here <p and ip are the scaled quantities. Since scaled quantities also satisfy the original equations, Eq. (|3.ti|l fixes 

e = 1 . (3.8) 

The scaling relations (|3.7|) along with (|3.8|l and successive transformations, first t — > 1 and then L — > t, gives scale 
invariant solution of the form 



h(t,z) = r a (j>(l,z/t) , 

f{t,z) = t-^(l,z/t) . (3.9) 

The equations above express arbitrary solution in terms of initial data (at t = 1). As stated earlier we work with 
t = 1 as our initial time and evolution is in sense of scaled time Lt with L > 1. 

From simple dimensional analysis we would have (3 = and a = ±1/2 (introducing a dimensional constant 
multiplying either dx or dy). However, as we will see later this would lead to a trivial solution: the homogeneous 
Kasner metric. The fact that we recover Kasner solution as a fixed point of a general inhomogeneous G2 metric though 
is non-trivial and it elucidates the appearance of this solution in earlier studies; it is used as a seed metric in different 
solution generating techniques, and more importantly, is known to describe "generic" cosmological singularity in the 
analysis of Belinskii et al. 22]. In order to get a non-trivial structure in fix-point analysis we need to analyze the 
anomalous dimensions for both the functions / and h. It is well known [2^| that the anomalous dimensions a =^ ±1/2 
and /3 7^ are fixed by initial or boundary conditions. In our case, cosmological vacuum solutions, initial condition 
refers to geometry at t = 0. Since, metrics of the type l|3.4|) are singular at t = there is not strict functional 
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constraint on the system described by l|3.5[l and (|3.6|l . Nevertheless, interpreting a posteriori the fixed points we can 
give some hints to understand the role of anomalous dimensions a and (3. 
Denoting L = exp(r) the RG equations are 

>(*)(£). ( , 10 , 

£-H)* + * + *(S)' + iW- 

We now investigate the fix-point structure of the equation set above, i.e., 

T = ^ \-0=(£)'+(£)' + U£fX - (3-12) 
rfr 4 I V*/ \<f> J % \ i> $ J 

System decouples to give 



-\/^%T' ( 3 - 13 ) 



(3.14) 



where A = a 2 + (3 — 1/4. The real solutions correspond to A > for z 2 > 1 and A < for z 2 < 1. The equations 
(|3.14|) and l|3.14|l can be easily integrated to give the fixed-points 



4>* = {z + ^z 2 - 1) ± ^ S , (3.15) 
^* = c,(z + V^ 2 - l)^ 2a ^ ( z 2 - , (3.16) 

were Cf is an integration constant. We have dropped constant of integration from ip* since this can be absorbed 
simply by scaling of x and y. Also, we note here that the new parameter A relates our spacetime metric with Kasner 
metric in the following way. The Kasner metric can be written in the form 

ds 2 = t^-^idz 2 - dt 2 ) + t l+d dx 2 + t^dy 2 , (3.17) 

where parameter d can be chosen positive or negative. A direct comparison with the metric (13.411 using (|3.9() gives 
the Kasner relationship 

4a 2 + 4/?-l = 0. (3.18) 
Therefore, A = case corresponds to the Kasner models. When A > (z 2 > t 2 ) the critical solution is: 

h(t,z) = r^z + Vz 2 -; 2 )^, 

/(*, z) = C/t («±V^+l/2)(a±v^-l/2) (z + v /^— ^)T2a^ ^2 _ f 2yA (3.19) 

and when A < (z 2 < t 2 ) the critical solution is: 

i z 

h(t, z) — t~ a exp(=Fy — A arccos -) , 

z . 



f(t,z) = c' f t a + A -V4^2 _ z 2)-A exp ( ±2ctV CA arccos -) . (3.20) 

t 

The spacetime is split into two regions separated by the light-cone z = t. In each region the solution takes one 
of the forms given by the expressions above. For A = we have the Kasner metric. In the general case (A ^ 0), 
solution (|3.19|l is the soliton metric that has been obtained earlier by the inverse scattering transformation with real 
poles from the Kasner metric, and the solution (|3.2U|) is the cosoliton solution generated also from the Kasner metric. 
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Both solutions have been studied in j24j (and references therein). Depending of the parameters, the light-cone z = t 
is singular for the solution (|3.19|) but is always singular for l|3.2U[l . This means that even though the metric and 
its first derivatives are continuous across the light-cone, the solutions cannot be matched across the light-cone (for 
a discussion on the matching of these metrics see [24]]). We would like to note here that the ± sign in the above 
equations is actually related to the symmetry under x <-> y, and as we will see later in the stability analysis same 
results hold for both the signs. Moreover, we would like to stress here that the solutions which we have obtained are 
actually the future asymptotic states, i.e., to which spacetime "prefers" to settle down. This makes this analysis very 
powerful since not only we recover in a very simple fashion a whole class of scale invariant solutions but also, due to 
"universality" , these solutions actually are the preferred asymptotic states. 

We will consider now the linear stability analysis. Since these solutions emerge as fixed-points of the RG map, the 
stability of these solutions is the stability of these fix-points. Let us define 



(3.21) 
(3.22) 



where 8<j> <C 1 and Sip <C 1. The above form of perturbations is chosen to facilitate the analysis since both unperturbed 
functions <fi* and tp* can diverge at z = ±1 and z — ±oo. 
The perturbation equations takes the form 



(h 
dSip 
~dr 



a + 2z- 



dz 2 dz 



<t>* 



2 | a + A + 2az— } — - + a— — 
dz dz 



We shall compute normal modes assuming the following form for the perturbations: 

54> = e WT p{z) , Sip = e UT a(z) , 
where w is a constant. The perturbation equations (|3.23(l can be written as: 



a + 2z- 



\°' = 



-up , 



2 1a + A + 2az—— I p — aa = ——uo 



From these equations it is easy to see that: 



uig = — 2y^Av z 2 — lp' — 2aojp 
Taking the derivative of this equation and substituting in (|3.25|) we get: 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



(1 - 2u)fJ + -p = 



z 2 - 1 



z 2 - 1' 



(3.27) 



Thus, the problem of solving linear perturbations around the fixed-points has been reduced to finding solutions of the 
above equation, and l|3.26|l completes the solution. The general solution of l|3.27|l is given by: 



p(z) = (z 2 - ir/ 2 +V4 (c lP »_+Jf {z) + c 2 Q^\z)) , 



(3.28) 



where Pjf and are Legendre functions of first and second kind, respectively, and c\ and C2 are arbitrary constants. 
From (|3.26() the solution for a can be easily obtained: 



(3.29) 



(z) = - 2v /A(z 2 - iyl*-V* [ Cl (zP^ 2 (z) P^ 2 1/2 {z)) + c 2 (zQ"X)i\z) - Q"l 2 1/2 (z)) 
(z 2 - ir /2+1/4 (ciP^/ 2 /2 (z) + c 2 Q^ 2 (z)) . 



-2a( 



We require that the perturbations be regular for all z. Since the differential equation (|3.27|l has four regular singular 
points at z = ±l,±oo we must consider three different regions: 1 < z < oo, — oo < z < —1 and —1 < z < 1. The 
first and the second region are equivalent. Let's first analyze the behavior in the first region. It is easy to see that 
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the leading terms of the two linearly independent solutions p\ and p 2 of l|3.27[l in the neighborhood of z = oo are 
respectively p\ « z u + 0(z ul ~ 2 ) and p 2 ~ px l nz + 0(z"~ 2 ). Both solutions are regular if a; < 0. Furthermore, from 
(13.26(1 . a is also regular provided that u> < 0. At z = 1 the leading terms of the two linearly independent solutions 
for p are: (z - l)"+ 1 /2 + q((z - 1)"+ 3 / 2 ) and constant + 0((z - 1)). For er the leading terms are: (z ~ l) u and 
constant + 0((z — 1)). Since uj is negative the first independent solution must be neglected so that we have regular 
solutions. For instance, tacking C2 = in l(3.28|l the solution is regular. 

Making the scale transformation given by l|3.7|l - H3.9|l we get the the behavior of the solution near the fixed-point : 

s^ = r P ( z /t) , 

Sifj = t u a(z/t) , (3.30) 

with u> < and z > t. Since the solution above is valid in the region z > t, taking limit t — > oo means z — ► oo 
as well. We can therefore distinguish two different situations. First, when we approach infinity in such a way that 
we encounter the light-cone z = t without crossing it. In this case asymptotically z/t = 1 and, bearing in mind the 
behavior of the regular solutions of (|3.27(l close to the point z = 1, the above perturbations behave as: 

S<j) ss Sip « constant x t u . (3.31) 

In the second situation we approach infinity through a path that does not encounter the light-cone. In this case 
z/t tends to a constant greater than 1, and the perturbations behave in a similar way that that given by the above 
expressions. So, in both cases the perturbations tend to zero when t — > oo which means that the fixed solution outside 
the light-cone is stable. 

Let's perform this analysis in the region — 1 < z < 1. The leading terms of the independent solutions of 13.27[) close 
to z = 1 are those described earlier for the same point changing z — 1 by 1 — z. When z = — 1 the behavior is the 
same changing z — 1 by z+ 1. From this we can not obtain a definite conclusion about the values of uj which make the 
solution regular. To do that is convenient to transform the equation l|3.27|l to a hypergeometric differential equation 
with coefficients a = b = 1/2 and c = 1/2 — uj. It is not difficult to see that one of the independent solutions the 
solution of H3.27fl and (|3.26|) is given by: 

pW = Cl(1+2) .Wi,i ; i_„ ; ;Lii 



a{z) = -20x^3^1 - z 2 (l + zy- 1 ' 2 

UJ 



1, /111 1- z\ 1 , N /3 3 3 l-z 
{bJ+ 2 )F 2'2 ; 2-^— -4-8^ (1 + z)F 2'2 ; 2-^^r 



-2aci(l + z)" +1/2 F^ 



111 l-z\ 



The form of the second independent solution depends on the values of u>. In any case, the above solution is regular 
when to is positive, except when u) = n+1/2 (n a positive integer). Therefore, with \z\ < t, the modes with uj positive 
dominate and the solution in the interior of the light-cone is not stable. 



IV. DISCUSSION 



In this paper we have shown how the RG method can be used to obtain the asymptotic regime of cosmological 
solutions. We have first illustrated the method applying it to homogenous cosmological models with a scalar field. 
These models have been extensively analyzed using, mainly, the so called "expansion normalized variables" ; which are 
a set of dimensionless variables . In those works the role of self-similar solutions in describing the asymptotic regime 
has been stressed. Since a fixed-point of the RG transformation is a scale-invariant solution, it is reasonable to think 
that the RG method will naturally render the same results (and we have manifestly shown this here, section [nj as 
those obtained using the "expansion normalized variables" . It is worthy to stress also that those normalized variables 
are well adapted to the dimensional and similarity analysis . These results bring the asymptotic behavior of the 
homogeneous cosmological models in a new perspective. Moreover, its simplicity as well as systematic approach show 
that it is a practical idea to implement RG for studying asymptotic behavior of spacetimes in General Relativity. 

We have also applied the RG method to diagonal, vacuum inhomogeneous G2 metric. In this case, the system 
reduces to a set of coupled partial differential equations whose analysis using the RG method is not difficult. The 
fixed-point is an exact solution depending on two parameters. This solution belongs to the class of soliton solutions. 
Soliton solutions are intended as those solutions which can be obtained by the inverse scattering transformation from 
a known one |20| . In particular, the found fixed-point is a soliton solution with real poles with origin at z = 0, whose 
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"seed" metric is the Kasner metric (soliton origin marks the origin of the light-cone z 2 = t 2 ). Since the metric i|3.4|) 
is invariant under a z translation the fact that the origin is at z — is not important. Moreover, a more general class 
of soliton solutions is that corresponding to a sum of solitons each with a different origin. Since the RG method gives 
the long-time behavior of the solution, the difference in the origins tend to zero as t tends to infinity leaving only 
one origin. It is interesting to note that, generically, the solution doesn't homogenize, let alone isotropize, for the 
final state is inhomogeneous. There is, however, a few particular cases for which the metric tends to a homogeneous 
solution: when A = the solution is the Bianchi type I Kasner metric, and for a particular value of the parameters a 
and A the solution is Ellis and McCallum family of vacuum Bianchi models 25]. Finally we would like to note that 
the fixed-points recovered in the RG technique gives all the scale-invariant exact solutions of vacuum G2 cosmologies. 

Now we can understand the relation of the anomalous dimensions a and (3 with initial conditions: exponent a gives 
the Kasner parameter of the "seed" metric to generate the soliton solution by means of inverse scattering technique. 
The parameter A (which fixes (3) determines the number of solitons in our solution. So, the two "initial conditions" 
used in the inverse scattering technique, number of solitons and the Kasner seed metric, fixes the anomalous dimension 
of system comprising of equations l|3.5J) and QH.fty . 

The RG method allows us to study the linear perturbations around the fixed-points as well. We have shown that the 
solution out-side the light-cone, which corresponds to the soliton solution with real poles, is stable against bounded 
perturbations, contrary to the solution inside the light-cone (cosoliton solution) which we find to be unstable. Let's 
note that although the stability analysis doesn't have any apparent dependence on the A parameter one should be 
very cautious to extend these results to the case A = 0. This case corresponds to the Kasner solution and this is the 
only case for which the fixed-point results in <f>*' = tp*' = 0, that makes the system of equations for the perturbations 
singular. 

There remain a few open questions that we hope to deal- with in future works. First is related with the role of 
soliton solutions with complex poles. These can be obtained as a complex translation along the z axis of the soliton 
solution with real poles and they are regular in all the spacetime, even on the light-cones. Since this complexification 
doesn't alter the structure of the field equations it would be worthy to investigate whether the complex pole soliton 
solutions can represent an asymptotic state of the inhomogeneous metrics. Second, it would be interesting to extend 
the analysis performed in this paper to general inhomogeneous non-vacuum metrics. Solutions with a scalar field are 
of particular interest in order to investigate issues like isotropization, scaling solutions, etc.. Finally this technique 
can be used to analyze long time behavior of spacetimes with more that four dimensions, for example, brane-world 
scenarios. 
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